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Abstract 

Graph pebbling is a network model for studying whether or not a given supply 
of discrete pebbles can satisfy a given demand via pebbling moves. A pebbling move 
across an edge of a graph takes two pebbles from one endpoint and places one pebble 
at the other endpoint; the other pebble is lost in transit as a toll. It has been shown 
that deciding whether a supply can meet a demand on a graph is NP-complete. The 
pebbling number of a graph is the smallest t such that every supply of t pebbles can 
satisfy every demand of one pebble. Deciding if the pebbling number is at most k is 
ilg-complete. 

In this paper we develop a tool, called the Weight Function Lemma, for computing 
upper bounds and sometimes exact values for pebbling numbers with the assistance 
of linear optimization. With this tool we are able to calculate the pebbling numbers 
of much larger graphs than in previous algorithms, and much more quickly as well. 
We also obtain results for many families of graphs, in many cases by hand, with much 
simpler and remarkably shorter proofs than given in previously existing arguments 
(certificates typically of size at most the number of vertices times the maximum degree), 
especially for highly symmetric graphs. 

Here we apply the Weight Function Lemma to several specific graphs, including 
the Petersen, Lemke, 4 th weak Bruhat, Lemke squared, and two random graphs, as 
well as to a number of infinite families of graphs, such as trees, cycles, graph powers of 
cycles, cubes, and some generalized Petersen and Coxeter graphs. This partly answers 
a question of Pachter, et al., by computing the pebbling exponent of cycles to within an 
asymptotically small range. It is conceivable that this method yields an approximation 
algorithm for graph pebbling. 
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*Much of this work was carried out while vising the Centre de Recerca Matematica during the research 
program on Variational Analysis and Optimization, Fall 2010. 
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1 Introduction 



Graph pebbling is like a number of network models, including network flow, transportation, 
and supply chain, in that one must move some commodity from a set of sources to a set 
of sinks optimally according to certain constraints. Network flow constraints restrict flow 
along edges and conserve flow through vertices, and the goal is to maximize the amount 
of commodity reaching the sinks. The transportation model includes per unit costs along 
edges and aims to minimize the total cost of shipments that satisfy the source supplies and 
sink demands. At its simplest, the supply chain model ignores transportation costs while 
seeking to satisfy demands with minimum inventory. The graph pebbling model introduced 
by Chung j7] also tries to meet demands with minimum inventory, but constrains movement 
across an edge by the loss of the commodity itself, much like an oil tanker using up the fuel 
it transports, not unlike heat or other energy dissipating during transfer. 

Specifically, a configuration C of pebbles on the vertices of a connected graph G is a 
function C : V(G)— >N (the nonnegative integers), so that C(v) counts the number of pebbles 
placed on the vertex v. We write \C\ for the size C(v) of C; i.e. the number of pebbles in 
the configuration. A pebbling step from a vertex u to one of its neighbors v reduces C(u) by 
two and increases C(v) by one (so that one can think of it as moving one pebble at the cost 
of another as toll). Given two configurations C and D we say that C is D- solvable if some 
sequence of pebbling steps converts C to D. In this paper we study the traditional case in 
which the target distribution consists of a single pebble at some root vertex r (one can peruse 
[TH IT5| ITT] for a wide array of variations on this theme). We are concerned with determining 
tt(G, r), the minimum number t of pebbles so that every configuration of size t is r-solvable. 
Then the pebbling number of G equals n(G) = max r n(G, r). Alternatively, n(G) is one more 
than the maximum s such that there is some root r and some size s configuration C so that 
C does not solve r. The primary focus of this paper is to exploit this duality with newly 
discovered algebraic constraints. 

1.1 Calculating Pebbling Numbers 

Given a graph G, configuration C, and root r, one can ask how difficult it is to determine if C 
solves r. In [18] it was determined that this problem is NP-hard. Subsequently, [191 [22] proved 
that the problem is NP-complete, with [19] showing further that answering the question "is 
ir(G) < fc?" is IT^-complete (and hence both NP-hard and coNP-hard, and therefore in 
neither NP nor coNP unless NP = coNP). Finding classes of graphs on which we can answer 
more quickly is therefore relevent, and there is some evidence that one can be successful in 
this direction. Besides what we share in this introduction, we show later that many graphs 
can have very short certificates that 7t(G) < k. 

The r-unsolvable configuration with one pebble on every vertex other than the root r 
shows that tt(G) > n, where n = n(G) denotes the number of vertices of G. In |20j it is 
proved that graphs of diameter two satisfy 7t(G) < n + 1, with a characterization separating 
the two classes (Class means n(G) = n and Class 1 means ir(G) = n + 1) given in [5j [8]. 
One of the consequences of this is that 3-connected diameter two graphs are Class 0. As an 
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extension it is proved in [TT] that 2 2d+3 -connected diameter d graphs are also Class 0, and 
they use this result to show that almost every graph with significantly more than n(n\gn) l l d 
(for any fixed d) edges is Class 0. Consequently, it is a very (asymptotically) small collection 
of graphs that cause all the problems. 

Knowing the pebbling number of a graph and actually solving a particular configuration 
are two different things, as even a configuration that is known to be solvable (say, one of size 
equal to the pebbling number) can be difficult to solve. Evidence that most configurations 
are not so difficult, though, comes in the following form. The work of [2] shows that every 
infinite graph sequence Q = (G\, G2, ■ ■ ■ , Gn, ■ ■ •) has a pebbling threshold Tg : N— >N, which 
yields the property that almost every configuration Cn on Gn oi size |Cjv| 3> t(N) is solvable 
(and almost every configuration of size |Cjv| "C t(N) is not). In papers such as [3j [9j [10] 
we find that Tg(N) is significantly smaller than tt(Gn) — for example, y/N as opposed to 
iV for the complete graph Kn, and roughly N2 v ' lgN as opposed to 2 N ~ 1 for the path Pn- 
Moreover, the proof techniques reveal that almost all of these solvable configurations can be 
solved greedily, meaning that every pebbling step reduces the distance of the pebble to the 
root. So the hardness of the problem stems from a rare collection of configurations. 

With these results as backdrop, [T] presents a polynomial algorithm for determining the 
solvability of a configuration on diameter two graphs of connectivity some fixed k. Further- 
more, [21] contains an algorithm that calculates pebbling numbers, and is able to complete 
the task for every graph on at most 9 vertices. Also, the proof in [7J that the d-dimensional 
cube is of Class is a polynomial algorithm (actually bounded by its number of edges nig n). 
Along these lines, our main objective is to develop algorithmic tools that will in a reasonable 
amount of time yield good upper bounds on ir(G) for much larger graphs, and in particular 
decide in some cases whether or not a graph is of Class 0. 

This latter determination is motivated most by the following conjecture of Graham in [7J. 
For graphs G and H, let GnH denote the Cartesian product whose vertices are V(G°H) = 
V(G) x V(H), with edges (u,x) ~ (v,x) whenever u ~ v in G and (u,x) ~ (u, y) whenever 
x ~ y in H . 

Conjecture 1 (Graham) Every pair of graphs G and H satisfy 7r(GaH) < 7r(G)n(H) . 

The conjecture has been verified for many graphs; see [13j for the most recent work. 
However, as noted in [16] , there is good reason to suspect that LpL might be a counterexample 
to this conjecture, if one exists, where L is the Lemke graph of Figure [TJ Since L is Class 0, 
Graham's conjecture requires that LnL is also, but it is a formiddable challenge to compute 
the pebbling number of a graph on 64 vertices. One hopes that graph structure and symmetry 
will be of use, but purely graphical methods have failed to date. The methods of this paper 
represent the first strides toward the computational resolution of the $64 questions, "Is 
n(LnL) < 64?". Certainly, these methods alone will not suflicell, but if they produce a 
decent upper bound then the methods of [21] might be able to finish the job. 

1 Yes, I'll pay if you beat me to it! 

2 We obtain evidence that n(LaL) < 108 in Theorem [TOl — in fact, for one root r we show Tr(LaL, r) < 68. 
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Figure 1: The Lemke graph 

1.2 Results 

The main tool we develop is the Weight Function Lemma El This lemma allows us to define a 
(very large) integer linear optimization problem that yields an upper bound on the pebbling 
number. This has several important consequences, including the following. 

1. The pebbling numbers of reasonably small graphs often can be computed easily. More- 
over, it is frequently the case that the fractional relaxation suffices for the task, allowing 
the computation for somewhat larger graphs. 

2. It is also common that only a small portion of the constraints are required, expanding 
the pool of computable graphs even moreJl One can restrict the types of constraints 
to greedy, bounded depth, and so on, with great success, seemingly because of the 
comments above. Potentially, this allows one to begin to catalog special classes of 
graphs such as Class 0, (semi-)greedy, and tree-solvable. 

3. The dual solutions often yield very short certificates of the results, in most cases 
quadratic in the number of vertices, and usually at most the number of vertices times 
the degree of the root. These certificates are remarkably simple compared to the usual 
solvability arguments that chase pebbles all over the graph in a barrage of cases. One 
can sometimes find such certificates for infinite families of graphs by hand, without 
resorting to machine for more than the smallest one or two of its members. This was 
our approach in Section I3.3[ for example. 

4. Our method gives trivial proofs of 

(a) 7r(C 2fc ) = 2 k and n(C 2 k+i) = 2[2 k+1 /3\ + 1, which we write as \(2 k+2 - l)/3] , and 

(b) C„ is Class for k > n/2(lgn — lglgn), where G^ denotes the k th graph 
power of G (as opposed to the Cartesian power G k ). This answers a question of 
|20j . who defined the pebbling exponent of G minimum such e = e 7r (G) for which 

3 We present some findings along these lines in Section 15721 with graphs on 15 and 20 vertices. 
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^(G^) = n(G). Thus e 7r (G ! ) < n/2(lgn — lglgn) (see Theorem [T6l). which is 
fairly close to the obvious lower bound of n/21gn. 

In this paper we apply the Weight Function Lemma to several specific graphs, including 
the Petersen, Lemke, 4 th weak Bruhat, Lemke squared, and two random graphs, as well as to 
a number of infinite families of graphs, such as trees, cycles, graph powers of cycles, cubes, 
and some generalized Petersen and Coxeter graphs. 

2 The Weight Function Lemma 

Let P n be the path V\V 2 • • -v n on n vertices. Then n(P n ) = 2 n_1 is easily proved by induc- 
tion. In particular, any configuration of at least 2™" 1 pebbles solves v%. But one can say 
more about smaller v i-solvable configurations as well, with the use of a weight function w. 
Define w on V{G) by w{y n -i) = 2\ and extend the weight function to configurations by 
w(C) = J2vev w{v)C{v). Then a pebbling step can only preserve or decrease the weight of 
a configuration. Since the weight of a configuration with a pebble on v\ is at least 2 n_1 , we 
see that 2 n_1 is a lower bound on every t^-solvable configuration. In fact, induction shows 
that every v i-unsolvable configuration has weight at most 2 n_1 — 1, which equals Y27=2 w ( v i)- 
That is, this inequality characterizes fi-unsolvable configurations on P n . The Weight Func- 
tion Lemma [2] generalizes this result on trees, and we explore the applications of the lemma 
in the following sections. 

2.1 Linear Optimization 

Let G be a graph and T be a subtree of G rooted at vertex r, with at least two vertices. For a 
vertex v £ V(T) let v + denote the parent of v; i.e. the T-neighbor of v that is one step closer 
to r (we also say that v is a child of v + ). We call T a strategy when we associate with it a 
nonnegative, nonzero weight function w with the property that w(r) =0 and w(v + ) = 2w(v) 
for every other vertex that is not a neighbor of r (and w(v) = for vertices not in T). Let 
T be the configuration with T(r) = 0, T(i>) = 1 for all v £ V(T), and T(t>) = everywhere 
else. With this notation note that the path result above can be restated: C is fi-unsolvable 
if and only if C(vi) = and w(C) < w(T), where T is the strategy T = P n with associated 
weight function w. 

Lemma 2 [Weight Function Lemma] Let T be a strategy of G rooted at r, with associated 
weight function w . Suppose that C is an r-unsolvable configuration of pebbles on V(G). Then 
w{C) < tu(T). 

Proof. By contrapositive and induction. The base case is when T is a path, which is proved 
above. Suppose w(C) > w(T), let y be a leaf of T, and define P to be the path from y to r 
in T, with P y being the subpath from y to its closest vertex x on P of degree at least 3 in 
T (or r if none exists). Denote by T' the tree T — P y + x, and among all such r-unsolvable 
configurations, choose C to be the one having largest weight on T'. The restriction w' of w 
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to T" witnesses that T' is a strategy for the root r, so induction requires that w'(C) < w'(T'). 
Likewise, the restriction w y of w to P y witnesses that P y is a strategy for the root x. Because 
w(C) = w'(C) + w y {C) and w(T) = w'(T') + w y (P y ), we must have w y (C) > w y (P y ) which 
by induction means that the restriction of C to P y — x solves x (and so x 7^ r). Let C x be 
the resulting configuration after moving a pebble from P y — x to x. Since C is r-unsolvable, 
so is C x . Now w(C x ) = w(C), but w'(C x ) > w'(C), which contradicts the initial choice of 



For a graph G and root vertex r, let T be the set of all r-strategies in G, and denote by 
Zg,t the optimal value of the integer linear optimization problem Pg )r : 



Max. \,C(v) s.t. w(C) < w(T), and T G T with witnessing weight function w . (1) 



We also let Zq,t be the optimum of the relaxation, which allows configurations to be ra- 
tional. We will find the relation Za, r < L^c,rJ useful at times. The following corollary is 
straightforward. 

Corollary 3 Every graph G and root r satisfies n(G,r) < zq^ t + 1. 

Proof. By definition, the pebbling number is one more than the size of the largest unsolvable 
configuration. □ 

Until now, one could only use trees in an individual manner: n(G,r) < n(T,r) for 
every spanning tree T rooted at r. The Weight Function Lemma allows one to consider all 
subtrees rooted at r (not only spanning trees) simultaneously, which we will see is significantly 
more powerful. One strength of the method is that the relaxation frequently has an integer 
optimum. This means that the dual solution will point out which tree constraints certify the 
result, and because the dual problem has only n(G) — 1 constraints there are at most that 
many such trees in the certificate. Experience has shown, however, that usually one can find 
a certificate with only deg(r) trees (or sometimes a few extra). We will see this behavior 
starting in Section [21 

2.2 Basic Applications 

We begin with the pebbling number of trees, whose formula was first discovered and proved 
in [7]. View a tree T with root r as a directed graph with every edge directed toward r. 
Then a path partition V of T is a set of edge-disjoint directed paths whose union is T. One 
path partition majorizes another if its nonincreasing sequence of path lengths majorizes that 
of the other. A path partition is maximum if it majorizes all others. We can use Corollary 
[3] to give a new proof of the following result of [7J. 

Theorem 4 For a tree T and root r we have vr(T, r) = XIpgp ^ ep — l^l + l; where ep denotes 
the length (number of edges) of P. 



C. 



□ 
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Proof. We begin by showing that a maximum size r-unsolvable configuration has pebbles 
on leaves only, and in fact on all leaves. Indeed, if C has a pebble on the nonleaf x, then we 
define a pushback of C at x to be any configuration obtained by removing the C (x) pebbles 
from x, adding 2C(x) + 1 pebbles to one of the children of x, and adding 1 pebble to all other 
children of x. Certainly, if C is r-unsolvable and has no pebbles past x (on the subtree of T 
rooted at x, minus x itself), then the pushback will also be r-unsolvable, and thus satisfy the 
constraints of P^ r It will also be larger than C. The configuration C* that places 2 ep — 1 
pebbles on the leaf of the path P 6 V is one possible result of pushing back the empty 
configuration, and so satisfies the constraints of Pt,v Hence ir(T,r) > J2pev^ ep ~ l^l + 1- 
For the upper bound we prove that C* is optimal by using induction to show that the 
optimal configuration has 2 ep — 1 pebbles on the leaf yp of P for every P EV. This is true 
if \V\ = 1, so suppose \V\ > 2 and let Q denote one of the paths in V whose leaf has the 
highest weight in w, with a tie going to one of the shortest length. This is to guarantee 
that the graph T — Q + x, where x is the root of Q, is a tree (i.e. is connected). In order 
to maximize the number of pebbles that satisfy ^2 P£V w(yp)C(yp) < J2 P( z V w(P) we would 
transfer as many pebbles from yQ to other leaves as possible because their weights are at 
most w(yp) and we could add extra pebbles when the weight is smaller. But by induction 
on T — Q + x we know from the constraint Ylpev-{Q} w (.yp)C{yp) < J2per-{Q} W {P) that 
each C(yp) < 2 ep — 1, with equality for all P ^ Q if and only if C is maximum. Therefore, 
since w(P) = w(yp)(2 ep — 1) for all P G V, we have 

w(y P )C(y P ) < 5>(P)- J2 ^(?/p)(2 ep -l) 
Pev Pev-{Q} 

= w{Q) 

= w(y Q )(2^-l), 

which implies that C(yp) < 2 e< 3 — 1. □ 

A slight weakness of these tree constraints is that they do not classify unsolvable config- 
urations on trees the way that they do on paths. This is because they let in a few solvable 
configurations. For example, consider the star on four vertices with one of its leaves as root 
r. Then the configuration with 2 pebbles on each of the other two leaves is r-solvable and 
satisfies all tree contraints. Since it is the average of the two r-unsolvable configurations 
that place either 1 and 3 or 3 and 1 on those other leaves, it cannot be cut out by the tree 
constraints that don't cut out at least one of these two other constraints. In this case it 
doesn't hurt us, since the strategy bound yields n < 5 and the actual pebbling number is 
5, but it can cause trouble on graphs in general. For example, we know that the 3-cube Q 3 
in Figure [2] has pebbling number 8, so that the shown configuration C is solvable (pebbles 
from to top must be split in two directions in its solution). However, no strategy recognizes 
its solution, and Corollary [3] yields only n(Q 3 ) < 9 (the three rotations of the strategy in 
the center in Figure [2] certify this). One can see where the aforementioned star appears in 
the Figure [2] configuration on Q 3 and is exploited accordingly: moving pebbles from the 5 
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r r r 

Figure 2: A solvable configuration (left) not recognized by any tree strategy; a canonical 
strategy (center) used for certifying n(Q 3 ) < 9; a simplified nonbasic strategy (right) that 
does the same 

along one edge yields a (3, 1) configuration, while splitting the moves along two edges yields 
a (2,2) configuration. 

3 General Applications 

In this section we illustrate the method more fully by presenting short proofs of both known 
and new results. We begin by relaxing strategies in the following way. We now use the 
term basic to describe the strategies as currently defined. A nonbasic strategy will satisfy 
the inequality w(v + ) > 2w(v) in place of the equality used in a basic strategy (see Figure 
[2]). The following lemma shows that we can use nonbasic strategies in an upper bound 
certificate since they are conic combinations of a nested family of basic strategies. Thus the 
use of nonbasic strategies can simplify and shorten certificates significantly. 

Lemma 5 If T is a nonbasic strategy for the rooted graph (G,r), then there exists basic 
strategies T±, . . . , TJ. for (G, r) and nonnegative constants c%, . . . , Cfe so that T = Yli=i c (^i- 

Proof. We use induction, as the result is true when T has two vertices since T is basic 
then. Given T, let S be a basic strategy on the edge set of T, define c to be the largest 
constant for which cS < T, and denote T' = T — cS. Then some vertex v of G satisfies 
cS(v) = T(v), so V has fewer vertices than T. Also, because S is basic, any vertex u whose 
unique ur-path contains v also satisfies cS(u) = T(u), which means that V is connected, 
and hence a strategy. Moreover, T' is nonbasic since every nonneighbor x of r has T'(x + ) = 
T(x + ) — cS(x + ) > 2T(x) — 2cS(x) = 2T'(x). By induction, T' is a conic combination of 
basic strategies, and so therefore is T. □ 

We use conic combinations of strategies to derive, for some a, the inequality |C| = 
Y2v^=r C{v) < a for r-unsolvable configurations C . From this we surmise that 7r(G) < [ctj +1- 
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Figure 3: Petersen Class certificate 



Instead of writing our strategies algebraically, it will be somewhat easier to show them 
graphically. We will display them so as to derive fn^2 v ^ r C(v) < J2v^r m ^( v ) < ma f° r 
some sequence {m v } v with m = min^m^, and let the reader divide by m. In fact, in many 
instances we will derive m v = m for all v ^ r, which makes for the following observation. 

Lemma 6 [Uniform Covering Lemma] Let T be a set of strategies for the root r of the graph 
G. If there is some m such that, for each vertex v ^ r, we have ^2 Te q-T(v) = m, then 



3.1 Specific Graphs 

It has been said in jest that every graph theory paper should contain the Petersen graph, so 
we get it out of the way first. 

Theorem 7 Let P denote the Petersen graph. Then n(P) < 10. 

Of course, the vertex lower bound implies 7r(P) = 10, but since the focus of this paper 
regards upper bounds, we prove them only. 

Proof. The 3 strategies shown in Figure [3] certify the result. □ 

Without such nice symmetry, the Lemke graph requires a different certificate for each 
possible root. 

Theorem 8 Let L denote the Lemke graph and suppose r ^ v\. Then 7r(L,r) < 8. 
Proof. We show the strategies for each root vertex in turn, below. 



7r(G,r) 



n{G). 



□ 



For r = f 2 : 



9 



For r = f 3 : 




Without the edge i>4i>8, the resulting graph would view v 3 and v 4 symmetrically. Using 
that symmetry, the solutions for r — v 3 become solutions for t> 4 . 
For r = f 5 : 
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□ 



Because of the 3-cube-like configuration with (1,1,1,5) on (t> 5 , t> 6 , t> 7 , v$), the best that 
our tree strategies can muster is ir(L,vi) < 9. Thus, to show that L is Class one must 
handle v i by more traditional methods. 

To illustrate that larger graphs can be tackled, we move on to one of order 24 that 
is not Class 0. Define the (weak) Bruhat graph of order m (see Figure H]) to have all 
permutations of [m] as vertices, with an edge between pairs of permutations that differ by 
an adjacent transposition. One can recognize it as the Cayley graph of S m , generated by 
adjacent transpositions, and also note that B4 is the cubic Ramanujan (expander) graph of 
[6] . Intuitively, expander graphs would seem to have low pebbling numbers, but because B4 
has diameter 6 we have the lower bound ^(B^) > 64. We give here a fairly tight bound. 

Theorem 9 Let B4 be the Bruhat graph of order 4- Then T[{B/±) < 72. 

Proof. Because the graph is vertex transitive, only one root must be checked. The 3 strategies 
shown in Figure H] certify the result. We combined them into one figure, separated by edge 
styles. □ 

Next we consider LnL, the square of the Lemke graph L. Because L has diameter 3, LnL 
has diameter 6, and so n(LnL) > 64. Strategies deliver the following upper bounds. Since 
the bounds are not that tight, we do not pursue bounds for all 64 vertices, although it is 
likely (since v 1 is the most problematic root in L) that the upper bound of 108 works for all 
roots (7*1, r 2 ). 

Theorem 10 Let LnL be the square of the Lemke graph L. Then 



Proof. For r = v 1 one can verify that a quarter of the sum of the four strategies in Figure [5] 
yields the weights 




108 if r = vx , 
96 if r = v s , and 
68 if r = t> 4 . 
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Figure 4: Bruhat graph £> 4 (left) and strategies (right) 
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giving the bound of 63 + 4(7 + 3 + 1) + 1 = 108. One verifies that each Tj is a strategy by 
making sure that each nonzero entry has a corresponding entry in its column or row with at 
least twice its weight that is joined to it by an edge in the appropriate copy of L. 

We use a similar argument when r = v$, using the strategies from Figure El along with 
their transposes, and divide by 8 to obtain the bound 63 + 8(3 + 1) + 1 = 96. 

Likewise, for r = v 4, the strategies and their transposes from Figure [7] yield an upper 
bound of 63 + [8(.6)J + 1 = 68. □ 

One can see where the improvement from 108 to 96 comes from. We obtain a bound from 
a sum of strategies by dividing by the minimum weight, and since the graph has diameter 
six, the maximum weight is at least 32. So we need to increase the minimum weight as much 
as possible by including more strategies, but there are diminishing returns. In the case of 
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Figure 5: Strategies for LnL at (i>i,i>i). 
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Figure 6: Strategies for LnL at (v 8 ,v 8 ). 
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Figure 7: Strategies for LaL at (t> 4 ,t> 4 ). 
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root (vx,Vi), it has four neighbors, and so a fifth strategy increases the maximum weight to 
at least 64. This is why the root (ug, ug), having eight neighbors, fares better. On this basis 
one might expect the root (t>4, U4) to perform best. In fact, it does even better for a different 
reason: it has eccentricity 4. It is amusing that the order of this graph is out of the range of 
computing, while these strategies were found fairly easily by hand. 

3.2 Random Graphs 

Given our introductory comments and observations, we also tested some modest sized ran- 
dom graphs. Let R15 denote the 15-vertex graph given by the adjacency list 

([2, 4, 5, 6, 12, 13], [1, 3, 4, 8, 11, 12, 14], [2, 4, 6, 7], [1, 2, 3, 5, 7, 9, 14], [1, 4, 6, 8, 11, 15], 
[1, 3, 5, 9, 13, 14], [3, 4, 11, 15], [2, 5, 10, 13, 14, 15], [4, 6, 10, 11], [8, 9, 11], 
[2, 5, 7, 9, 10, 12, 15], [1, 2, 11, 13], [1, 6, 8, 12], [2, 4, 6, 8], [5, 7, 8, 11]). 

We generated this graph with edge probability .35. It has diameter three, with 10 of the 
vertices (including v$) having distance two from all others. We found that -R15 is at least 
nearly Class as follows. Our simple certificate of 3 basic strategies for the root r = v 9 is 
below. 

2xi + 2x 2 + 2x 3 + 4x 4 + 2x 5 + 2x 7 + x 8 + x 12 + x 13 + 2x u + x% 5 < 20 
2x 2 + 2x 7 + 2x 8 + Ax w + Ax u + 2x 12 + x l3 + 2x 15 < 19 
2xi + 2x 3 + 2x 5 + Ax e + x 8 + x 12 + 2xi 3 + 2x u + x 15 < 17 

We note that there are 20, 422 basic strategies of depth at most two rooted at r. These were 
generated by our java code in a few minutes on my Core 2 Duo 2.66GhZ PC, with 2GB 
RAM and 250GB HD running Ubuntu linux 9.04. CPLEX software solved the resulting 
linear optimization problem instantly, delivering a dual certificate of a conic combination of 
11 of the strategies. By hand, we grouped them according to their neighbor of r and were 
able to combine each into a single nonbasic strategy. By trading weights between some of 
the strategies, we arrived at the simplified certificate above. The certificate might be more 
easily viewed in the table of coefficients, below, in which the sum of the constraints appears 
below the line. Thus 7r(i? 15 ,t> 9 ) < 15. 
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The same upper bound of 15 was obtained for all other roots (see the Appendbo) except 
for viq, for which we found 7r(_Ri 5 , v 10) < 16. Note that -R15 contains the 4-cycle (t>4, v 3 , vq, vq), 

We should point out that whenever there are too many possible strategies to compute in a reasonable 
amount of time and space — on random graphs with 30 vertices and diameter 4 it was typical to crash 
memory after a day of running — one can generate a large collection of strategies at random rather than by 
exhaustion and obtain optimal results far more quickly. This was a useful speed up even for R20 with roots 
V3, v w , V12, vir, and v 2a . 
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with vg adjacent to the root v w - - the same dreaded configuration we discussed in Q 3 
obstructs us here: (1,5, 1,0). More to the point, let C be the size 15 configuration (1, 1, 5, 
1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1) on Ri 5 . Because v 3 is distance 3 from the root vi and C is 
empty on the neighbors of t>io, the only way to move a pebble to the root involves splitting 
the pebbles from v 3 as in the dreaded configuration solution. Hence no strategy can detect 
the Vio-solvability of C . However, we can blend strategies with some case analysis and slight 
amount of old fashioned analysis as follows. 

Theorem 11 The graph R± 5 is Class 0. 

Proof. As we have already shown that 7r(i?i5,r) = 15 for all r except for viq, it suffices 
to prove that 7r(i?i5, v iq). Let C be an r-unsolvable configuration, where r = v\q. We 
first consider the case that C contains a pebble on a neighbor of r. In this case we have 
x§ + xg + Xn > 1, and the CPLEX certificate 
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shows that |C| < 14. In this format, the left column holds the multipliers of each constraint, 
and the bottom row is the result of the linear combination of them. Division by 995 yields 
the result. 

Thus we may assume that x§ + Xg + xu = 0. Next we consider the case that at most 5 
pebbles are distance 3 from r; that is, x-y + x 3 < 5. Here we have the certificate 



/ 


84 


1 





1 








































168 























1 


1 





1 
















6 


2 


4 


2 


2 


3 


2 





8 











2 


4 


4 





33 




13 


2 





2 


4 





4 


1 





8 











2 


2 





25 




29 


2 





2 


4 


2 


4 


2 





8 











1 


2 





27 




55 


2 


4 


2 


2 





2 


2 


8 











2 


3 


4 


4 


35 




23 


2 


4 


2 


2 


4 


2 


1 


8 











2 


4 


4 


4 


39 




42 


2 





2 


2 


4 


2 


4 











8 


4 


2 





1 


31 




6 


2 


4 


2 


2 


4 


2 


4 











8 


4 


2 


2 


1 


37 




36 


2 


4 


2 


2 


4 


2 


3 











8 


4 


2 


2 


4 


39 


\ 


504 


504 


504 


504 


504 


504 


504 


840 


504 





840 


504 


504 


504 


504 


7476 / 



On the other hand, at least 7 pebbles at distance 3 from r yields the following certificate. 
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Hence we may also assume that X\ + £3 = 6, and so there are exactly 9 pebbles on the 9 
vertices at distance two from r if we assume that \C\ = 15. 

Now, v\ and V3 can both reach all neighbors of r, so if either x\ > 4 or 23 > 4 then we 
can place a pebble on any neighbor of r. This forces each vertex at distance two from r to 
have at most, and hence exactly 1 pebble. Thus the splitting of the 4 pebbles will place a 
pebble on r. This contradiction means that X\ = £3 = 3. 

Similarly, then, the pair of v 1 and t>3 can place a pebble on any neighbor of r via a 
common neighbor, again forcing 1 pebble on every vertex at distance two from r. However, 
this allows v\ and V3 to each place a pebble on the same neighbor of r via disjoint paths. 

This final contradiction shows \C\ < 14 again, finishing the proof. □ 



Next we considered the random graph i? 20 having adjacency list 

([6, 8, 11, 12, 14, 15, 16, 17], [4, 5, 6, 7, 8, 10, 15, 16, 17, 18, 19, 20], 

[4, 6, 8, 12, 14, 20], [2, 3, 5, 6, 8, 9, 12, 15, 18, 19], [2, 4, 7, 12, 14, 15, 16, 18, 20], 

[1, 2, 3, 4, 7, 8, 14, 15, 19], [2, 5, 6, 8, 11, 12, 13, 14, 15, 17, 18], 

[1, 2, 3, 4, 6, 7, 10, 11, 14, 15, 17], [4, 10, 11, 13, 14, 17, 19, 20], 

[2, 8, 9, 16, 18, 19, 20], [1, 7, 8, 9, 13, 14, 16, 18, 20], [1, 3, 4, 5, 7, 13, 16], 

[7,9,11,12,19,20], [1,3,5,6,7,8,9,11,18], [1,2,4,5,6,7,8,19], 

[1, 2, 5, 10, 11, 12, 18, 20], [1, 2, 7, 8, 9, 19], [2, 4, 5, 7, 10, 11, 14, 16, 20], 

[2, 4, 6, 9, 10, 13, 15, 17, 20], [2, 3, 5, 9, 10, 11, 13, 16, 18, 19]), 

generated with edge probability .4. The graph has diameter two, does not have the form 
of the Class 1 characterization, and so must be Class 0. Indeed, from the 9, 371 basic 
strategies of depth at most two for the root r = v\ (other roots had between 8,000-28,000 
such strategies), CPLEX delivered the following certificate involving 16 of them (which we 
won't bother to simplify) in Figure El 

Theorem 12 The graph R 2 o is Class 0. 

Proof. The certificate for ^(i^c^i) = 20 is shown in Figure [HJ The certificates for all 
remaining roots are shown in the Appendix. □ 



3.3 Graph Classes 

Next we turn our attention to classes of graphs, and begin with an extremely simple proof 
of the pebbling numbers of cycles, first proved in [20] . 

Theorem 13 For k > 1 we have 7r(C 2k ) = 2 fc and 7r(C 2 fc+i) = \{2 k+2 - l)/3] . 

Proof. For both results we use two basic strategies: one path in each direction. For even 
cycles the paths of length k will overlap in the vertex opposite the root. This yields 7r(C2k) < 
2(2 fc — l)/2 + 1 = 2 k . For odd cycles with k > 3 the paths will be of length k + 3, which 
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Figure 9: The generalized Petersen graph Ps i2 . 

gives n(C 2k+1 ) < L2(2 fc + 3 - l)/(2 2 + 2 3 )J + 1 = [(2 k+2 - l/2)/3j + 1 = \(2 k + 2 - l)/3]. For 
k < 2, paths of length k + 1 suffice. □ 

Now we consider a generalization of the Petersen graph. For m > 3 and d > 2 define 
P m; d to have vertices it and v^x, where 1 < i < m and X is a binary fc-tuple for < k < d. 
Furthermore, uv^ is an edge for all i, and v^xv^x- is an edge for all i and nonempty X, 
where X~ denotes the truncation of X obtained from dropping its final digit. Finally, for 
every i and length d — 1 X, we include the edge i^x^+i.x (addition modulo m), with the 
exception that when i = m — 1 we use f m _i,x^o,x+ instead, where X + denotes the (d — 1)- 
tuple that satisfies N(X + ) = N(X) + 1, and N(X) is the natural number represented by X 
in binary. Figure M shows the graph P 52 ; it is easy to check that P3 2 is the Petersen graph 
P. Also, P m ,d has diameter 2d when m > 3 and 2d — 1 when m — 3. 

Interest regarding these graphs comes from two sources. In [20] the problem is raised 
of finding the minimum number e(n) of edges in a Class graph on n vertices. Because 
of the Petersen and Wheel graphs, we have e(n) < 2n — 2. Blasiak, et al. [4J, show that 
e(n) > [_3rz / 2J , and conjecture that e(n) = 3n/2 + o(l). Furthermore, they believe that the 
graphs P m ,d are all Class which, if true, would prove the conjecture because it is 3-regular 
except for the central vertex u, having degree m (with n = m [2 d - 1) + 1). The graphs also 
appear in [12] as n-vertex graphs having the fewest edges among those of minimum degree 
3 and radius d. In fact, the intuition for the conjecture comes from this result. 

Here, we are most interested in the graphs P m)2 for m > 3. Note that P m>2 has rotational 
symmetry: the exceptional edges for k = 2 can be "rotated onward" by swapping i> 10 with 
Vi t i in the drawing. This makes P m .2 transitive on the set {y^x} for fixed length X. Thus, 
when calculating vr(P m2 ), one only need consider the three vertices u, v = v\$, and w = i^o 
as root. 

Theorem 14 For all m > 4 we have the following, where n = n(P m 2 ) = 3m + 1: 
1- 7r(P mi2 ,it) = n, 
2. 7r(P mj2 , v) < n + 5, and 
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Figure 10: The graph Pg,2 and its main strategy T for root u. 




Figure 11: The four strategies for root v in Pg^. 



3. ir(P m> 2,w) <n+17. 

Proof. When the root is it, we use all rotations of the following basic strategy T (see Figure 
ITOj) . Let t>! have weight 4, t> 10 and V\ t x have weight 2, and each of their two neighbors have 
weight 1. Clearly, of the sum of all rotations of T has weight 4 everywhere but u, and so the 
Uniform Covering Lemma applies. 

When the root is v , we build slightly more complex strategies. First we use the nonbasic 
strategy S, having weight 3 on t> and v 0) i and weight 1 on each of their neighbors. Next, 
write the rotations of T as To, ... , T m _i. For j G {0, 1, 2}, we build the basic strategies Sj 
from these by combining all T (i ^ 0) with i = j mod 3, with the exception that if j = 1 
and % = m — 1 then we do not include the vertices vq$ and t>o,i twice. Moreover, each Sj 
includes weight 8 on u (see Figure flTl) . One quarter of the sum of these four strategies 
has weight on v, weight 6 on u, and weight 1 everywhere else. Now the Weight Function 
Lemma (actually Corollary [3]) applies. 

The description for the strategies when the root is w is almost identical to that for root 
v (see Figure [T2|) . In this case, one quarter of the sum of the four strategies has weight on 
w, weight 13.5 on v, weight 6 on u, and weight 1 elsewhere. □ 

Next we define the generalized Coxeter p-graph C(p) for odd primes p = 2q + 1. Like the 
graphs P m ,d, C(p) has the potential (for p > 5) to be a Class graph with fewer than In — 2 
edges. Set V { = | < j < p} and V = U q i=0 Vi, so that C(p) has n = p[q + l)= ( p ^) 

vertices. Define the 2pq = 2n — 2p edges + 1) and (0,j)(i,j) for each 1 < i < q and 

< j < p. The vertices in Vq have degree q and all others have degree 3, making C(7), the 
original Coxeter graph, 3-regular. Here we show the following n + 0(y/n) bound. 



20 




Figure 12: The four strategies for root w in Pg^- 



Theorem 15 For every prime p > 11 we have ir(C(p)) < n + 14g + 14. 

Proof. First we exhibit the symmetry of G = C(p). To do so, we note that the arithmetic 
that follows will be modular in q in the first coordinate, with the speciality that we use 
the representative q in place of 0, and modular in p (as normal) in the second. Let p be 
the automorphism that sends vertices (i, j) to + 1) for all < i < q and < j < p. 
This is a rotation of the vertices within each V^. Let a generate Z* and define a to be the 
automorphism that sends (0,j) to (0,aj) and to {i + l,aj) for all 1 < i < q and 

< j < p. While permuting V , the key aspect of a is that it also permutes the sets of 
vertices from V% to Vi + \. Together, p and a act transitively on Vq and on Vq. This means that 
we only need calculate n(G, (0,0)) and vr(G, (1,0)). (It turns out that C(7) is fully vertex 
transitive and so the root (0, 0) suffices.) 

Next we define the q basic strategies T, (1 < i < q) for the root r = (0,0), with each 
root child having weight 4. The child of r in Tj has weight 4 and two children, (i, i) and 
(i, —i). We describe the descendents of only, as those of (i, —i) are identical but with 
their second coordinate negated. The left subtree under is the path (i,2i), (i,3i), . . ., 
(i,qi)- The right subtree under (i, i) is a collection of q — 2 paths, starting with children 
(k,i) for 1 < k < q with k (jL {i, (p±i)/2}. The child (k,i) = (k,sk), for s = ik" 1 , and from 
it hangs the path (k, (s + l)k), . . ., (k, qk). 

Now we look at the sum of the weights of vertices over all strategies, and by the symme- 
tries described above we need only consider the vertices (i, j) for i 6 {0, 1} and < j < q. 
We note first that (0, j) is the grandchild of (j, 0) in Tj, and so has weight 1. The left path 
under (1, 1) shows that vertex (1, j) has weight 2 2_ - ? in Tj. For j > 3 we see that if k < j 
then vertex (1, j) has weight 2 fc ~- ? ~ 1 . Therefore, except for the q weight 4 children of r and 
their 2q weight 2 children, every other nonroot vertex has weight 1. This gives total weight 
(n — 1) + (4 — l)q + (2 — l)2g, implying that tt(G, r) < n + 5q. 

Consider the other root r = (1,0). Here we need many more basic strategies — we 
will define T« for i e {±1} U {2, . . . , q] — and the weights of the children of r equal to 
8 instead of 4. Similar to above, the strategy T_i will be identical to T« except that the 
second coordinates will be negated. We will describe the first four levels of each strategy 
explicitly, and their remaining levels implicitly. The child of r in Ti is (1, 1), having children 
(1,2) and (0, 1). From (1,2) hangs the path (1,3) and (1,4), while the children of (0, 1) are 
(2, 1), . . . (q, 1). The children of each (i, 1) are (i, 1 ± i), subject to them not having already 
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listed in some Tj, which we now describe. The child of r in Tj is (0, 0), which has the single 
child (z,0). The two children of (i, 0) are the two children of which are (z,±2z) and 

(0,±i), correspondingly, with the exception that (i, ±2i) does not appear in T g because it 
already appears in T-ti (as (g, ±1)). At this point, there are q + 1 vertices of weight 8, q + 3 
vertices of weight 4, and 4q — 2 vertices of weight 2, with all other vertices having weight 1. 

The remaining construction of the strategies proceeds in stages, with each new vertex 
added with several fractional weights adding up to 1. In the first stage, each vertex (a, b) 
that is adjacent to a current leaf (a, c) of some strategy is added to that strategy as a child 
of (a,c), accounting for weight 1/2. The other 1/2 weight comes from adding it as a child 
of (0, b) as well. In each subsequent stage, every vertex (a, b) that is adjacent to a current 
leaf (a, c) is added as a child of (a, c) in every strategy that contains (a,c), accounting for 
weight 1/2. The other 1/2 weight comes from adding it as a child of (0, b) as well. 

Hence we obtain ir(G,r) < (n - 1) + 7(q + 1) + 3(q + 3) + l(4q - 2) + 1 = n + Uq + 14. □ 

We note that the number of strategies for the root (1,0) can be reduced whenever two 
strategies (not including T±i) share no vertices (other than (0,0)). Thus one can define the 
intersection graph G p of the family of sets Wi, 2 < i < q, where Wi = V(Ti) — {(0,0)}. 
With chromatic number Xp — x(G p ), we obtain the upper bound of n + 7% + 7q + 21. For 
example, Xp 1 when p < 7. In those cases, further improvements can be made as well, 
and it is not too hard to show that 7r(C(5)) < n + 6 and ir(C(7)) < n + 15. 

Finally, we discuss powers of cycles. For a given graph G and integer k, we denote by G^ 
the graph on the same vertex set as G, with edges uv whenever the distance distc(w, v) <k 
in G. For example, G 1 -™ -1 ^ = K n for every connected G, where n = n{G). Pachter, et al [20J, 
define the pebbling exponent of G to be the minimum e = (G) for which G^ is Class 0. 
Consequently e 7r (G) < n — 1 for all G. The problem raised in [20] is to find e 7r (C n ). Here we 
prove the following. 

Theorem 16 The pebbling exponent of the cycle C n satisfies 

'2/1 < e , (Cn) < . 

Ig n Ig n — Ig Ig n 

Proof. The lower bound follows from the general fact that ir(G) > 2 diam ( G ) for all G, along 
with the observation that diam(Cn 6 ' ) ) = |~n/2e~|. Therefore, a requirement for Class is that 

For the upper bound, we prove that ir(Cn ) = n for n = (2k + l)2 fe + 3 = f(k). Then 
we show that ix(Cn ) — n for f{k — 1) < n < f(k). Our method will be to split the cycle 
into two identical paths with endpoints at the root r, and use identical strategies on each 
one. We will invoke the Uniform Covering Lemma M to obtain the result. 

In more detail, let us give a useful labeling of the vertices V of C n as follows. We partition 
V = uf+ 2 (C/i U Wi) so that 

• each Ui and Wi induces a path in C n , 
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• U = Wo = {r}, 

• Ui+i follows Ui when traversing C n clockwise from r, 

• Wi+i follows Wi when traversing C n counterclockwise from r, and 

• U k+2 = W k+2 . 

The two identical paths mentioned above are seen to be U = \Jg$Ui and W = U k ^W h 
where now we see that the term 'split' was a white lie because of the overlap on U k +2 = 
Wk+2- We will describe a family T of strategies on U with the property that, for some m 
(actually 2 k+1 ) and all u E Ufi^C/j, we have ^ TeT T(u) = m, while for all u E U k + 2 we have 
^ Tgr T(-u) = m/2. Then we copy these strategies symmetrically onto W and, because of 
the overlap, the Uniform Covering Lemma |6] applies. 

Next we describe the vertices within each U{. First, \U\ = 2 k for i E {l,k + 2}, and 
\Ui\ = 2 k - 2 k - i+1 for 2 < % < k + 1 (thus n = (2k + l)2 k + 3). We clockwise order the 
vertices of Ui = {vi >0 , . . . , Vi^-i} in their natural order by subscript, and will find it useful 
to identify Vij with the encoding [i, bj], where bj is the fc-bit binary representation of j (so 
that leading zeros are not suppressed). For example, v^q is encoded as [3,00110] in C355 
since k = 5 in that case. Also, the root r is encoded as [0, k ], where we write x J to denote 
the concatenation xx- ■ -x of length j. Furthermore, we use the notation Vj to mean some 
binary word of length j, subject to context. That is, [1, v^] is always a vertex in U\ but [2, v^] 
is not always in U2 — every vertex in Ui looks like [2, Ov^-i]. Moreover, for 2 < % < k + 1, 
vertices in Ui look like everything but [i, V^Vk-i+i]- 

Now we describe a partially ordered set P that we will use to define our strategies. The 
elements of P are the vertices U, and the covering relations are given by 

1. [0,0 k ] > [l,v fc ] for all v fc , 

2. [l,Vfc_ix] > [2, 0vfc_i] for all Vfc_i and all x E {0, 1}, 

3. [i, Vj_!Vfc_ja;] > [i + 1, v$_i&Vfc_$] for all 2 < i < k, all Vi_i 7^ l l_1 , all ^k-i-, and each 
x,6G{0,l}, 

4. [i, V- 2 0v k -ix] > [i + 1, l i_1 0vjfe_i] for all 2 < i < k, all v fc _i, and all x E {0, 1}, 

5. [k + 1, vjfc] > [A; + 2, v fc ] for all v fe 7^ l k , and 

6. [A; -h 1, 1*- 1 0] > [A; + 2,l d ]. 

All other relations of P are determined by transitivity. 

For z E P, define its downset D(z) = {y E P \ y < z). Notice that each D = D([l, v^]) 
forms a tree in P. Indeed, simple induction shows that if Vjt = V\ • ■ • Vk then every vertex 
of Ui+i in D has the form [i + l,Xjt>i • • -Vk-i], where Xj is anything but 1*. This means 
that [i + l,Xi ■ ■ ■ x(V\ ■ ■ ■ Vk-i] has exactly one element from D that covers it, namely either 
[i, x x --- Xi_iVi ■ ■ ■ v k _iV k ^i+i] or [i, x 2 --- x(o x ■ ■ ■ t> fe _;t> fe _ m ]. So no cycles exist in D. 
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For each v fc , then, define the basic strategy T = T(v fc ) to have vertices {r} U D([l, v fc ]), 
with edge zy whenever z > y is a covering relation in P. In order that T is a strategy in 

(2^ ) h 

C n we must verify that the distance d = d\st(z,y) between z and y in C n is at most 2 . 
When 2; = r, we have 2/ = t>ij for some < j < 2 k — 1, so d = j + 1. When 2 = vij, we 
have y = t>2,Lj72j, so d = 2 fc — |~j/2~|. When z = v^j and 2 < i < k + 1, then we can write 
j = h2 k ~ l+1 + / for some < / < 2 k ~ l+l . From relation [3] or H] we have y = t>2,j', where 
j' = h2 k ~ l+1 + b2 k ~ t + \)/2\ and b G {0, 1, 2}. Thus the greatest distance comes from relation 
SI in which case d = \U\ - j + f = 2 k - \l/2] . 

Note that the characterization of elements in D([l, v*.]) implies that each vertex of Ui is 
in 2 l_1 basic strategies when 1 < i < k + 1, and in 2 fc when i = fc + 2. Because each T is 
basic, this means that Y1t^( v ) = ^ fc+1 for all v G Ui (resp. Wi), 1 < % < k + 1, and 2 fc 
for t> G (resp. ^+2). Now the overlap from Uk+2 = Wk+2 gives sum m = 2 k+l for all 

v 7^ r, and the Uniform Covering Lemma [6] applies. 

Finally, whenever n is smaller than f(k) we simply erase sufficiently many vertices of 
Uk+2 = Wfc+2 (but don't renumber any indicies/encodings), as they are leaves in all strategies 
and so don't destroy the uniform covering. When such vertices are exhausted continue 
erasing vertices of Uk+i U Wk+i with similar results. Since f(k) — f(k — 1) = 2 fc+1 , no more 
considerations are necessary. □ 



4 Remarks 

In this paper we have shown several different strengths of the Weight Function Lemma 
in combination with linear optimization, highlighting its versatility. It has been used to 
compute upper bounds on and exact values of the pebbling number of small graphs. It has 
also been successful in calculating the pebbling numbers of much larger graphs than previous 
algorithms. For such graphs having too many strategies than time allows to construct, the 
technique of creating a smaller set of them at random seems to perform just as well. This 
is most likely due to the property that nonoptimal solutions (derived from having fewer 
constraints) seem, in most instances, to be near optimal (have the same floor function). In 
fact, by restricting strategies to be breadth-first search, one obtains upper bounds on the 
greedy pebbling numbers of graphs (which requires pebbling steps to move toward the root). 
The method also yields results for many families of graphs, in many cases by hand, with 
much simpler and remarkably shorter proofs than given in previously existing arguments. 
This is especially so with highly symmetric graphs. We note also that the technique can be 
used in conjunction with more traditional arguments, as in Theorem [TT], and it has delivered 
an array of upper bounds, such as n, n + c, and n + o(n), most of which are the best 
known and might possibly be best possible. It's two main shortcomings are the inability 
to overcome the kind of splitting structure found in cubes, for example, in which solutions 
to some configurations require nontree solutions, and the difficulty in dealing with large 
diameter, although success has been found with cycles and their graph powers, in addition 
to Petersen and Coxeter generalizations. When the technique gives upper bounds, it would 
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be of great use to know how good the bound might be. That is, does the Weight Function 
Lemma yield an approximation algorithm for graph pebbling? 

Question 17 Is there a constant c such that, for all graphs G and every root r G V{G), 
z G ,r<cn(G,r)? 

For example, is c = 2? The cube Q d shows that it couldn't be any smaller. 

Theorem 18 For all d > 1 and all r e V(Q d ) we have z Q d^ r < 2n(Q d ,r). 

Proof. We exploit the symmetry of Q d as follows. Because Q d is vertex transitive we need 
only consider one root r. We identify V(Q d ) with the power set of {1, ... , d} and take r = 0. 
We define a single strategy T for r and then apply every permutation of {1, . . . , d} to obtain 
other strategies. Finally we average over this collection of strategies. The result will be that 
n(Q d ) will be at most one more than the sum of the weights in this average. 

For each 1 < k < d such that 2 k < Q we define a k = [Q/2 k \ and b k = Q mod 2 k . 
We form T by first taking a neighbor of r. This neighbor will be a set of size 1, to which we 
assign the weight 2 d ~ 1 — this is step i — 1, with k = d — i. For future steps i, while 2 k > (f), 
we continue adding a single neighbor, a set of size i, to the current leaf of T, and assign the 
weight 2 k . If 2 k < (f\, however, we add a k + 1 vertices to the current leaves: a k of these 
will have weight 2 k and one will have weight b k . All of these will be connected to leaves of 
weight 2 k+1 and none will be connected to leaves of weight b k+1 . This is possible because 
the degree of each vertex of size i — 1 is k + 1, and a k < {k + l)a k+ i whenever k + 1 < d — 1, 
so there are enough potential neighbors to accomplish this. 

Hence we obtain weight 1 on average for vertices of size % for which 2 k < (^) , and weight 
27© otherwise. This yields the bound z Q ^ r < 1 + £to max{ g) , 2 k } < ELo (S + 
Yl~ = \2 k <2 d+l = 2ii{Q d ). □ 

If one restricts their attention to only polynomially many strategies, this linear optimiza- 
tion technique becomes a polynomial algrorithm. It would be useful to investigate how good 
the approximation can be under these circumstances. 
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6 Appendix 

6.1 Certificates for Ri$ 

We list the missing certificates in order, with the all- zeros column signifying the root. The 
format is the same as that for i? 15 at t> 10 in Section 13.21 No attempt was made to find the 
simplest set of strategies. 



577 







8 


4 















2 




4 







2 














2 


1 




2 


25 




40 







8 















2 







3 







2 




4 












4 




2 


25 




1250 










4 




1 


i 










4 









4 


2 

















3 




2 


27 




542 




















8 





2 




4 







2 




4 




1 




2 


2 




2 


27 




190 




















8 


4 







4 







2 




1 

















4 


23 




217 




















8 


4 


1 




4 







2 




4 









2 







4 


29 




62 










4 












8 


2 









4 


2 




1 




















21 




646 























8 







2 




4 


2 




1 









4 


4 







25 




94 







4 


2 















2 




2 




1 


2 




4 




8 







2 




2 


29 




91 







8 


4 












1 







4 




2 


2 




4 




2 




2 


4 




2 


35 




300 










2 









8 


4 


2 




1 




2 


2 




4 




2 




2 


2 




4 


35 




344 







4 















2 


2 




2 







2 




4 




8 




4 


1 




2 


31 




175 







4 


2 












1 












2 


2 




4 




8 




4 


2 




2 


31 




472 







4 


2 












1 


2 




2 




2 


2 




4 




8 




4 







2 


33 









10004 


10002 


10000 




9992 


9998 


9999 


10000 


10002 


10000 


9998 


10004 


10002 10001 


9998 


14000 


-J 












/ 


4 




























1 











1 


























4 


































1 





1 


























2 




4 








2 2 








(1 





1 


2 








11 




























2 










4 


2 


1 























7 




























1 










4 











2 











1 





7 




























1 










4 





2 





2 











1 





9 




























2 



















4 





2 





2 


1 


2 


11 




























2 













2 


2 





2 


2 


4 1 








2 


15 


























V 






8 





8 8 


8 8 


8 


8 


8 


8 


8 8 


8 


8 


8 


112 




















( 


2 





8 












2 







4 







2 


1 


4 




2 


C 





23 
















6 


2 


8 




















4 







2 


1 


4 




2 


4 


2 


29 














34 


2 


8 




















4 







2 


4 


4 




2 


3 


2 


31 
















23 


4 












8 


4 










4 




2 





1 




2 


4 





29 
















1 





















8 







4 




2 


1 


2 




4 


C 





21 
















42 





















8 




2 







2 


4 


1 







C 


4 


21 
















16 


4 












8 


4 







2 


4 




2 


1 


2 




2 


J 





33 
















3 


4 












8 


4 










4 




2 





2 




2 


3 


2 


31 
















19 


















4 


8 




2 


4 




2 





1 




4 


C 


2 


27 
















7 


4 















4 


8 




2 


4 




2 


1 


2 




4 


c 


2 


33 
















15 


4 















4 


8 







4 




2 





2 




4 


3 


2 


33 














\ 




336 336 









336 


336 336 336 




336 


336 




336 336 


336 


336 336 


336 


4704 


J 









27 



135 












2 













4 













1 













2 




9 \ 


31 





















4 










8 


4 




3 





2 












21 


2 









8 
















2 





1 


2 




4 


4 


2 







2 




25 


1 




2 




4 


8 








4 










1 







2 


2 


2 












25 


38 




4 













8 







2 


3 





2 




4 


2 


2 







4 




31 


10 




2 
















4 










8 


4 




4 


2 





1 




2 




27 


1 












2 








4 










8 


4 




4 


2 


2 


2 




1 




29 


103 









4 


2 





1 










4 





2 







2 


2 


8 




2 




27 


39 




8 













4 


4 





















2 


4 


1 




2 




25 


30 




8 




2 


1 








4 







2 













4 


4 












25 


49 




2 




8 


4 













2 


4 





2 







4 


2 







1 




29 


4 




2 




8 


4 








1 




2 








2 




4 


4 





4 









31 


54 












2 





8 


4 




2 


1 





2 




4 


2 


2 


2 




4 




33 


82 




2 




1 


2 








4 




2 


2 


8 


4 




4 


2 


2 












33 




996 


998 


1000 





995 


996 




998 


1000 


995 


996 




994 


996 


998 


999 


996 




13957 / 








10 












































1 




1 












6 


4 
































2 


1 










7 












2 





2 


2 


4 








1 




























9 














1 








2 


4 








2 




















1 







9 














3 








2 


4 








2 





2 














1 







11 














6 








2 








4 








1 











2 


2 







11 














2 























4 





2 








1 










7 














2 





2 

















4 





2 











2 


1 




11 














2 





2 

















4 





2 








2 


2 


1 




13 














6 





2 














2 





2 


2 


4 


2 








1 




15 














2 




















1 























2 




3 












V 




24 


24 


24 


24 





24 


24 


24 


24 


24 


24 


24 


24 


24 


24 


336 


J 






7 




8 










4 










2 





1 










4 


4 












23 


5 




8 













4 
























4 


3 









2 


21 


79 


























2 





8 


2 




4 


2 












1 


19 


1 









4 


2 


3 



























2 







8 







19 


14 




8 







2 


4 










1 













2 


4 


4 




2 







27 


120 




2 




4 


8 


3 










4 
















2 












2 


25 


31 




4 




2 








8 







2 





1 


2 




4 





2 









4 


29 


91 















4 


8 







2 


1 





2 




4 


2 












4 


27 


4 









2 





3 













2 


8 


4 




4 


2 












2 


27 


20 




4 
























4 





2 







1 


8 












19 


85 




4 
























4 





2 







2 


8 









1 


21 


104 









4 



















4 





2 







1 







8 







19 


37 















4 










2 


1 


8 


4 




4 










2 




2 


27 


7 




2 




4 


1 


4 










2 


4 





2 







2 







8 




2 


31 




1006 


998 


997 


999 


996 





9 


98 


1000 


998 


998 


996 


1002 


1001 


9£ 


8 


c 


!)- 




13985 



\ 



( 


99 


2 


4 


8 








4 





1 











2 


2 


2 





25 > 


26 


2 


4 


8 








4 





2 


2 








2 


2 


1 





27 




2 


4 








8 


4 


1 








4 


2 





2 


2 


4 





31 




59 























1 


2 


2 


4 


2 


1 








12 




51 


2 


4 




















3 


4 


8 


4 


2 


2 





29 




31 














4 








1 


4 


4 


8 


4 


2 








27 




12 


2 











4 








4 





2 








2 


1 


8 


23 




65 














4 


2 





4 





2 








2 


1 


8 


23 




48 


2 


1 








4 


2 





4 





2 








2 


2 


8 


27 




123 


4 


2 





8 


3 


2 





2 


4 


2 





2 


2 


4 





35 




13 


2 














2 





1 


4 


4 


8 


4 


2 








27 


V 


998 


998 


1000 


1000 


1001 


1000 





1000 


999 


998 


996 


998 


999 


999 


1000 


13986 j 



203 


























2 


4 


1 














7 


208 


2 














1 

















2 


4 








9 


7 








2 


4 


8 











2 





3 














19 


66 











4 





4 


1 





2 














8 





19 


37 


2 


3 


2 


4 





4 


2 




















8 





25 


116 








2 


1 








4 























8 


15 


38 


4 


8 


4 








2 


2 





2 


1 


4 


2 


2 


4 





35 


30 


4 





2 





8 


4 








2 








2 





1 





23 


5 


2 





2 


4 


8 


4 


2 

















1 








23 


77 


2 





2 


4 


8 





2 





2 


1 


4 


2 











27 


6 








2 





8 


4 


2 





2 





4 











1 


23 


37 





8 


4 








2 


2 





2 


2 


4 


4 


2 





1 


31 


6 


2 


8 


4 








1 


2 





2 





4 


4 


2 








29 


31 


2 


8 


4 


4 





2 


2 





2 





4 


4 








1 


33 




1000 


1007 


1004 


1008 


1000 


1002 


1004 





1002 


1001 


1004 


1002 


999 


1006 


1002 


14041 



28 



82 


2 





2 


4 








8 





1 














2 





19 x 
21 


15 








2 


4 








8 





2 














2 


3 


98 

















2 





1 











8 


4 








15 


3 























4 





2 











1 


8 


15 


12 























4 





2 








1 


2 


8 


17 


4 





8 


4 


4 





1 








2 








4 





4 





27 


38 


4 


8 


4 











1 


4 











4 


2 


4 


2 


33 


11 





8 


4 


4 





2 





4 


2 








3 





4 





31 


67 


2 





2 





8 


4 





1 


2 











2 


2 


4 


27 


58 


4 





2 





8 


4 


1 





2 














2 


4 


27 


14 


2 





4 


4 





1 


8 




















2 





21 


120 











2 





1 





4 


4 


8 








2 





2 


23 


64 


4 


8 


4 


4 





2 





1 


2 











2 


4 





31 


8 


4 


8 


4 








2 


2 


4 





1 





2 


2 


4 


2 


35 




998 


1000 


1000 


1000 


1000 


1000 


1000 


997 


1000 


998 





1001 


998 


999 


997 


13988 / 





28 


8 





2 


4 


3 





2 





2 
















2 





23 






4 


8 








4 


2 


3 


2 





2 
















2 





23 




3 














4 





4 





1 


4 


8 
















21 






4 








2 





3 





4 








4 


8 
















21 






27 


8 





2 


4 


4 


2 


2 





1 
















2 


2 


27 






33 





8 


4 


4 


1 


2 


2 


4 





2 













4 


2 


33 






52 














2 


1 


4 





4 


4 


8 













4 


27 






59 








2 





1 


4 





4 


2 


2 










8 


2 


2 


27 






26 





8 


4 


4 


2 


2 


2 


4 


2 


2 













4 


1 


35 




I 




472 


472 


472 


472 


472 


472 


472 


472 


472 


472 


472 







472 


472 


472 


6608 


) 


/ 


8 


8 





2 


4 








2 





1 
















2 





19 




10 




















2 





1 


2 


4 




8 








2 


19 




15 


8 





2 


4 


4 





2 






















2 


1 


23 






14 


8 








4 


4 





2 





2 
















2 


1 


23 






5 





2 


4 





3 


8 


2 





4 


2 


2 










4 





31 






13 





4 





2 


1 





2 


8 





4 


2 










4 


4 


31 






27 





4 


2 


1 








2 





2 


2 


4 




8 





2 


2 


29 






2 





2 


4 





4 


8 


2 





4 


2 


2 













1 


29 






17 








4 


2 


4 


8 


2 





4 


2 


2 













1 


29 






13 





2 


4 


1 


4 


8 


2 





4 


2 


2 













2 


31 






24 





4 


2 


2 


1 





2 


8 


2 


4 


2 










4 


4 


35 




V 


296 


296 


296 


296 


296 


296 


296 


296 


296 


296 


296 


2' 


fi 





296 


296 


4144 J 



19 


2 


4 


2 





4 








8 





4 


2 


2 








3 


31 




109 





8 


4 











2 








1 


4 


4 


2 





2 


27 


3 








2 


8 


3 





4 





4 


2 


2 














25 




22 














4 


1 


2 


8 








2 


2 


4 





4 


27 




27 


2 











1 








8 





4 








4 





4 


23 




30 














2 


1 





8 





4 





2 


4 





4 


25 




16 




















2 


8 


2 


4 








3 





4 


23 




3 


2 


8 


4 





2 


1 


2 


4 





2 





4 











29 




7 


4 





4 


8 


4 





1 


2 


4 


2 


2 











2 


33 




116 


4 








8 


4 





4 





4 


2 


2 


2 


1 





2 


33 




54 


4 





4 





2 


8 


2 





4 


2 


2 


2 


1 








31 




28 


2 





4 








8 


1 





4 


2 


2 


2 


4 








29 




35 


4 





4 





4 


8 


2 





4 


2 


2 





3 





2 


35 




8 





4 


2 








1 


2 


8 





4 





2 


4 





4 


31 


-) 




1002 


1004 


1004 


1008 


1006 


999 


1005 


1002 


1004 


1001 


1004 


1002 


1001 





1003 


14045 





( 


14 


4 


2 


2 


4 


8 


2 








1 








2 


2 








27 \ 




6 


4 


2 


2 


1 


8 


4 








2 








2 


2 








27 




1 


4 


2 





2 


8 


4 








1 








2 


2 


2 





27 




5 


4 





2 


2 


8 


4 








1 








2 


2 


2 





27 




2 


2 


1 


4 


4 





2 


8 


























21 




14 


2 


2 


4 


4 





1 


8 




















2 





23 




10 


2 





4 


4 





1 


8 





2 














2 





23 




2 


2 


4 


1 


2 





2 





8 


2 


4 





2 


4 








31 




24 


2 


1 











2 





8 


2 


4 





2 


4 


4 





29 




26 





4 


2 


1 





2 








4 


4 


8 


4 


2 


2 





33 


V 


208 


208 


208 


208 


208 


208 


208 


208 


208 


208 


208 


208 


208 


208 





2912 J 



6.2 Certificates for R20 

We list the missing certificates in order, as above. 
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